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We perform a detailed quantum dynamical study of non-equilibrium Josephson oscillations be- 
tween interacting Bose-Einstein condensates confined in a finite-size double-well trap. We find that 
the Josephson junction can sustain multiple undamped Josephson oscillations up to a characteristic 
time scale Tc without quasipartcles being excited in the system. This may explain recent related 
experiments. Beyond a characteristic time scale Tc the dynamics of the junction is governed by fast, 
quasiparticle-assisted Josephson tunneling as well as Rabi oscillations between the discrete quasi- 
particle levels. We predict that an initially self-trapped BEC state will be destroyed by these fast 
dynamics. 



One of the striking manifestations of quantum mechan- 
ics on a macroscopic level is the particle current induced 
by the phase difference betvireen two coherent wave func- 
tions connected by a weak link, known as the Joseph- 
son effect [3]. For tunneling between the macroscopic 
wave functions of two Bose-Einstein condensates (BEC) 
trapped in a double- well potential this phenomenon leads 
to temporal oscillations of the population imbalance z 
between the two condensates even in the ground state 
[2, 3]. However, unHke in superconducting Josephson 
junctions, the interaction between the atoms in the con- 
densates gives rise to regimes of fundamentally new dy- 
namical behavior. If the initial population imbalance 
exceeds a critical value depending on the interaction 
strength, the large-amplitude Josephson oscillations (de- 
localized regime) cease, and the BEC is trapped in one 
of the two wells with only small-amplitude oscillations of 
z about the non-zero mean value (self-trapped regime). 
This complex, non-linear dynamics has been theoretically 
predicted for the ground state 0, 0, @| and experimen- 
tally observed recently It is not only interesting in 
its own right but also relevant for any merging process of 
BECs, e.g., for quantum engineering or for producing a 
continuous source of condensed atoms [^, Q . 

In the recent experiment the Josephson junction is 
prepared in a non-equilibrium situation by ramping up 
the barrier between the condensates suddenly, in a non- 
adiabatic way. For this case an immediate damping of the 
Josephson oscillations has been predicted due to quasi- 
particle excitations which is, however, not observed 
in the small traps of Ref. revealing an incomplete 
understanding of the non-equilibrium dynamics. 

In this Letter we present a detailed study of the tem- 
poral non-equilibrium dynamics of Josephson-coupled 
BECs after non-adiabatically switching-on the Joseph- 
son coupling, including interatomic interactions as well 
as quasiparticle (QP) excitations. As the main result we 
find that in small traps multiple, undamped Josephson 
oscillations are possibe up to a time scale Tc. At this time 




FIG. 1: (Color online). A BEC in a double-well potential 
after abrupt decrease of the barrier height. The definitions of 
the parameters are explained in the text, see Eqs. ©-(O. 



scale the dynamics switches abruptly but continuously 
from slow Josephson to fast Rabi oscillations between 
the discrete QP levels. We also predict that the self- 
trapped behavior is destroyed by the Rabi oscillations, 
i.e. for time t > Tc the system switches to delocaHzed be- 
havior, if it has previously been in a self-trapped state. 
This highly non-linear behavior results essentially from 
a separation of energy scales in small traps with discrete 
QP level spacing A, which can be chosen larger than J. 
Switching on J lowers the ground state energy by the 
amount AE = J^iVi(0)A^2(0), where 7Vi(0), iV2(0) are 
the occupation numbers of the two BECs in the initial 
state at time T = 0. Thus, after a sudden switching 
two initially separated BECs are in an excited state AE 
above the coupled ground state. Because of the large val- 
ues of Ni , N2 this energy is sufficient to excite QPs out of 
the BECs. The time-dependent BEC amplitude acts as a 
perturbation on the QP system. However, transitions to 
QP states are not allowed in perturbation theory, because 
the frequency of the oscillations is less than their excita- 
tion energy, J < A. Our detailed calculations show that 
such transitions are only possible as a highly non-linear 
process after the characteristic time Tc. 

We consider a Bose-Einstein condensed atomic gas in a 



double-well trap as represented by Fig. [TJ Such a system 
is most generally described by the Hamiltonian 



where ^(r, t) is a bosonic field operator, and we assume a 
contact interaction between the bosons with g = Anas/m 
(fls is the s-wave scattering length). Vext is the exter- 
nal double-well trapping potential. Initially, the barrier 
between the two wells is assumed to be infinitely large, 
so that Josephson tunneling is absent. All bosons are 
condensed, and both condensates are in the equilibrium 
state. At time t = the barrier is suddenly lowered so 
that a Josephson weak link is established between the 
wells. This nonadiabatic process drives the system out 
of thermodynamic equilibrium. 

In order to develop the general non-equilibrium theory 
for this system and to analyze its dynamics, we wish to 
represent the Hamiltonian ([1]) in the complete basis of 
the exact single-particle eigenstates of the double-well 
potential Vext{^,t > 0) after switching on the coupling 
J. In this basis the field operator reads, 

^(r,t) = 0i(r)ai(O+02(r)a2(i) + ^¥'n(r)&„(i), (2) 

where (j>i (r) , (/)2 (r) are the respective ground state wave- 
functions of the two wells after lowering the barrier, and 
the aa are the corresponding, time-dependent condensate 
amplitudes (c- numbers), aa{t) — \/A^ae**"'*', a — 1,2. 
This semiclassical treatment of the BECs neglects phase 
fluctuations. It is justifled for the experiments where 
the BECs are initially produced with fixed phase rela- 
tion and the particle number is sufficiently large. The 
applicability of the semiclassical ap pro ximation has been 
discussed in detail in Refs. fiol. [T3| and has been 

tested experimentally in Ref. [1J|. The quasiparticle dy- 
namics will be treated fully quantum mechanically. The 
index n ^ enumerates the exact single-particle excited 
states, with (p„(r) and 6„(i) the corresponding eigen- 
functions and bosonic destruction operators, respectively. 
Note that by including excited states we go beyond the 
frequently used two-mode approximation [1, Q for the 
BECs. For simplicity we assume that the ground state 
energies of the two wells before mixing are equal, Eq = 0, 
and that the wavefunctions of the excited states extend 
over both wells. Inserting the field operator ([2]) into 
Eq. ^ and evaluating the overlap matrix elements in 
a straightforward way, the Hamiltonian takes for i > 



the form, H = Hbec 
condensate particles. 
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FIG. 2: (Color online). Time-evolution of condensed and non- 
condensed particles for the initial conditions z(0) = —0.6, 
^(0) = and interaction parameters m = u' = 5, j' = 60, 
k — Q (delocalized regime). 5 QP levels were included in the 
numerical evaluation, (a) The dynamics of the BEC popula- 
tion imbalance z(T') is shown (solid black line). The dashed, 
blue line shows, for comparison, the behavior without QP 
coupling, u = j' — k = 0, in agreement with Ref. (b) 
z(T) vs. 0{T) map. The arrow indicates the direction of time 
evolution. The time T = Tc is marked by the black dot. It is 
seen that at T = Tc the system changes dynamically from a 
S = 0toaS = 7r junction with more erratic phase evolution, 
(c) Time evolution of the non-condensed particle population. 
The black curve (also in the inset) shows the particle occupa- 
tion of the first level n[,^' , while the red curve is the sum of all 
five levels nt. For T < Tc the two curves practically coincide. 



with the Josephson coupling J and the interaction be- 
tween condensed particles, U > 0. Hqp corresponds to 
single-particle excitations. 



U' 



(4) 



where £"„ are the (bare) QP energies, and U' is the repul- 
sive interaction between non-condensed particles. Mixing 
between the BECs and the QP system is described by 
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(a*a2 + a'^aiYbl^bn + -(aja2&„6,i -|- h.c. 



(a>„)6j,6„ + -( 



(5) 
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^0 /. KO'a+— /(al^aaaaaa)- J {ala2+a2ai) , 



(3) 



Here the coupling constant J' arises as a QP-assisted 
Josephson tunneling as well as a pairwise QP cre- 
ation/destruction out of both BECs simultaneously. K 
represents the density-density interaction of condensed 
and non-condensed particles and the pairwise QP cre- 
ation/destruction out of each of the BECs separately. In 
deriving Eq. |[5]) we neglected the off-diagonal in n and 
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m elements because of different spatial dependence of the 
wavefunctions. 

To treat the non-equilibrium quantum dynamics of the 
system, we use the Keldysh Green's function G + C, gen- 
eralized to Bose-condensed systems. The QP part reads, 



{TMt)blit')) {TcK{t)hm{t')) 
{Tcbi{t)bHt')) {Tcbi{t)bm{t')) 



G„ 



Gn.' 



(6) 



where Tc impHes time ordering along the Keldysh con- 
tour, i.e. each of the normal and anomalous bosonic 
Green's functions G and F is a 2 x 2 matrix in Keldysh 
space. The condensate part is classical with trivial time 
ordering. 



CQ/5(t, t') 



aa{t)a*p{t') aa{t)ap{t') 
<{t)aW) <{t)ap{t') 



(7) 



The equations of motion for G -I- C are derived in a 
standard way Transforming the time variables to 

center-of mass and relative coordinates, T = {t + t')/2 
and T — {t — t'), observing that the Josephson dynamics 
depending on T is slow compared to the inverse QP ener- 
gies (J < A), the relative coordinate can be set t = in 
all propagators and self-energies. We treat the QP inter- 
action in Eq. (0]) within the self-consistent BogoHubov- 
Hartree-Fock approximation This will be sufficient 
for the present purpose, since QP collisions, neglected 
here, will play a role only for sufficiently high population 
of QP levels (see below). The normal and anomalous 
QP self-energies, I]„(r = 0,r), rj„(T = 0,T), are then 
diagonal in the QP level index n and read. 



S„(T) 



K{Ni + N2) + J'{ala2 + a;ai) 
2jC/'^G<„(T), 



(8) 
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with E = E and fl = il* . After lengthy but straightfor- 
ward calculations one arrives at the coupled set of equa- 
tions for the non-condensate propagators G'^{t = 0,T), 
F<(t = 0,T) and the complex condensate amplitudes 
ai(T), 02 (T), 

d_ 



i^G<m = r!„(r)F:„(r)-f]„(r)F<„(T), 
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2E„(T)j^^<„(T) 

f7„(r)G<„(r), 
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(10) 



^Qf^iiT) - [U\ai{T)\^ + KN.iT)] a^T) - Ja2{T) 



+J'Nb{T)a2{T) + i 



Eal{T)+'^al{T) 
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The equation for 02 (T) is obtained from Eq. (fTTI) by 
fli ^ 02- From Eqs. (fTO|) . (fTTI) '^^ compute the occupa- 
tion numbers for bosons out of condensate, Nj^^^T) — 
(6t(T+)6„(T)), NbiT) =. En^i'^hT), the condensate 



population imbalance z{T) — [Ni{T) — N2{T)]/N, nor- 
maHzed by the total particle number N — Ni{0) + 
-^2(0) + iV(,(0), and the time evolution of the phase differ- 
ence 0{T) = 02(T) — 9i{T). To absorb the large factors 
of particle numbers appearing in Eqs. l(9|)- lfTT|) it is use- 
ful to define the dimensionless parameters u = NU/J, 
u' = NU'/J, f = NJ'/J, k = NK/J, and n["^(r) = 
iVf")(T)/7V. 

Without coupling to the QP excitations {j' — k — Q), 
Eq. ifTTj) reduces to the two-mode model of Smerzi et al. 
[H|, exhibiting the self-trapped and delocalized regimes, 
with a Josephson oscillation frequency of 



UJ 



(0) _ 



= 2|J|v/l + w/2 



(12) 



in the linear regime (Eq. (10) in Ref. 01). When, how- 
ever, / 7^ 0, A: 7^ and QPs are excited, Nb{T) > 0, 
the QP-assisted Josephson tunnehng term in Hmix be- 
comes active (J'-term in Eq. (fTT]) ). One then expects an 
enhanced Josephson frequency, with roughly J replaced 
by J[l — j'ni,{T)] in Eq. (fT2l) . At the same time, Rabi 
oscillations of the nI^"'\t), i.e. of QP pairs between the 
BEG and the excited levels, with frequencies ujr ~ 2En 
set in, c.f. Eq. (fTOj) . As a result, in this QP-dominated 
regime one expects complex, high-frequency anharmonic 
oscillatory behavior. 

The complete numerical solutions of Eqs. (|9|)- lfTT|) for 
a finite-size trap with iV = 5 • 10^ particles (level spacing 
A = En+i — En = 10 J, taking 5 QP levels into account 
[l3|) are shown for typical parameter values in Fig. [2] 
for the delocalized regime and in Fig. [3] for the initially 
self-trapped regime. 




FIG. 3: (Color online). Same as in Fig. 2, but for u — u' — 25 
(self-trapped regime) and j' = 30. In (a) the dashed, blue line 
shows the behavior without QP interactions, u' = 0. 
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FIG. 4: (Color online). Dependence of the inverse time-scale 
on two main parameters j' and k, for the initial conditions 
as in Fig. 3. In the left panel the white area corresponds to j' 
and k values for which Tc was not found; we assume Tc ^ oo 
in this case. The right panel shows a collapse of all t^^ curves 
onto a single one with the simple law ~ J —k. The scatter 
is due to the ambiguity in the numerical definition of Tc. 



The results reproduce the expected behavior discussed 
above in the regime -with finite QP population. The pa- 
rameters were chosen such that the energy A_E stored 
in the two BECs by the initial, non-adiabatic switching- 
on of J is much greater than the QP energies, Ai? w 
Ja/7Vi(0)A^2(0) > A. The most striking and most im- 
portant feature seen in both figures is that nevertheless 
multiple undamped Josephson oscillations occur for an 
extended period of time without QPs being excited. The 
reason for this behavior is that in the initial state the 
QP population nb{t) is vanishing and, therefore, the QP- 
assisted Josephson tunneling term J' in Eq. ^ does not 
contribute. Hence, the Josephson oscillations have the 



bare frequency loj 



'■ to 



(0) 



< 2A which is not sufficient to 



excite a QP pair perturbatively. Only for times greater 
than a characteristic time Tc the highly non-linear dy- 
namics of the system makes QP excitations possible. In 
this long-time regime the finite QP population rn,{t) and 
fast oscillations of the BEC population imbalance z{t) 
stabilize each other mutually: nb{t) > impHes a QP- 
enhanced Josephson frequency, and the resulting fast os- 
cillations {ll>j > 2A) of z(t) can efficiently excite QPs via 
the mixing Hamiltonian |[5]). 

The fast, QP-induced dynamics impHes two further 
features. (1) As seen from Fig. [3] (a) an initially self- 
trapped state is destroyed and the system changes to a 
delocalized state at the same time when nb{T) becomes 
sizeable. (2) At the onset of the fast dynamics the system 
changes from a6' = Oto6' = 7r Josephson junction, see 
Figs. [2] (b), [3] (b). This can be understood qualitatively, 
in that the large phase difference 9{T) « tt is required to 
sustain the large Josephson current in the state with fast 
dynamics. 

Since the transition to the QP-dominated regime is not 
described by a Fermi golden rule, it is hard to analyse the 
time scale Tc analytically. We defined Tc numerically as 



the scale where ni,{t) first exceeds 0.05 and extracted it 
from our solutions. As seen in Fig. [3] (a), Tc is essentially 
independent of the QP interaction u' . The dependence 
of 1/tc on the parameters / and k (Fig.lH) is remarkably 
linear, and for j' < k no transition to the QP-dominated 
regime is found. 

To conclude, we have presented a detailed quantum 
dynamical study of the non-linear Josephson dynamics 
of BECs confined in a finite-size double-well potential, 
including coupling to quasiparticle states. Remarkably, 
the system can sustain multile, undamped Josephson os- 
cillations for an extended time period before quasipar- 
ticles get excited and the behavior changes abruptly to 
a regime of fast Josephson and Rabi oscillations. Only 
in this quasiparticle-dominated regime we expect strong 
damping of the oscillations due to inelastic quasiparticle 
collisions, equilibrating the system at a finite tempera- 
ture. This will be a subject of further research. The 
sharp but continuous transition from the Josephson- to 
the quasiparticle-dominated regime should be experimen- 
tally observable. 
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